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ABSTRACT 

Preceding page blank not filmed 

The conventional solutions for increasing the flutter 
speeds of particular configurations have included, struc- 
tural modifications, mass balancing, and geometry changes, 
which generally result in increased structural weight and 
reduced aerodynamic efficiency. Flutter suppression by ac- 
tive control of aeroelastic modes with closed loop feedback 
of appropriate parameters is a relatively recent approach 
which can significantly reduce these design penalties. This 
thesis is concerned with the development of organized synthe- 
sis techniques, using concepts of modern control theory, for 
the design of active flutter suppression systems for two and 
three-dimensional lifting surfaces, utilizing a control mo- 
ment gyro (CMG) to generate the required control torques. 

Incompressible flow theory is assumed, with the un- 
steady aerodynamic forces and moments for arbitrary airfoil 
motion obtained by using the convolution integral based on 
Wagner’s indicial lift function. Wagner's function is ap- 
proximated in the analysis by a sum of exponential terms and 
two-dimensional strip theory is used to model the finite wing. 
This formulation enables Laplace transform techniques to be 
applied to the system of equations. Additional states, based 
on second and third order time derivatives of bending and 
torsional displacements, can be identified and thus allow 
the system to be placed into standard state variable form, 
with the state transition and control distribution matrices 
being real. 

Using pole-assignment (arbritrary dynamics) techniques, 
the control system characteristics of the two-dimensional 
(typical section) wing are first examined to establish basic 
design trends with the results presented as a function of 
velocity on root locus diagrams. Linear optimal control 
theory is applied to find particular optimal sets of gain 
values which minimize a quadratic performance function. 

The closed loop system’s response to impulsive gust distur- 
bances and the resulting control power requirements are in- 
vestigated, and the system eigenvalues necessary to mini- 
mize the maximum value of control power are determined. Full 
state feedback with scheduled gains is required to obtain 
this minimum value. Reduced feedback designs with constant 
gain values, which do not significantly increase power re- 
quirements, are also determined. 

The vibrational properties of the wing in developing 
the three-dimensional analysis are described by a superposi- 
tion of the ’’uncoupled" modes of pure bending and torsion. With 
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full state variable feedback, complete modal controllability 
(ie.j arbitrary relocation of the system eigenvalues) of 
the system can be obtained with only a single CMG controller 
regardless of the number of assumed modes. The significant 
feedback parameters needed to obtain this control are deter- 
mined as a function of structural and important CMG para- 
meters and velocity. As was learned in the typical section 
analysis, the location of all the closed loop eigenvalues 
proved critical to particular designs, not only in meeting 
power level requirements, but also in influencing the occur- 
rence of higher mode instability. 

Conditions for system observability are discussed for 
both the typical section and three-dimensional wing cases. 
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I. INTRODUCTION 


1 . 1 Background 

Flutter can be defined as a self -excited, oscillatory 
instability of an elastic body in an airstream. Physically, 
this instability can occur when energy from the airstream 
is absorbed by the body more rapidly than it is dissipated 
by damping. The onset of the phenomenon is typically char- 
acterized by a flutter frequency and a critical airstream 
speed. Flight at airspeeds in excess of this critical flut- 
ter speed can result in catastrophic structural failure. 

The problem of increasing the flutter speed of particular 
configurations has in the past been generally ( solved at the 
expense of increased structural weight, through structural 
modification and mass balancing, or by reducing aerodynamic 
efficiency, through geometry changes. Weight penalties in 
the range of 10 to 25,000 pounds, for example, have been 
estimated (1,2) to meet flutter requirements in SST class 
vehicles. 

In more recent years, a substantial amount of experience 
has been obtained in developing active control systems for 
gust alleviation and the suppression of low frequency, in- 
sufficiently damped structural modes of an aircraft. Refs. 

3 and 4 are excellent survey papers of the work done in these 
areas. This experience, coupled with the potential perfor- 
mance payoffs available (1-5), has made the active suppres- 
sion of the unstable flutter modes both feasible and desir- 
able. Much of the work done in this latter area has occurred 
in the last 5 years, and some recent studies can be found in 
Refs, 6 through 11. In Ref. 6 a study, including analog sim- 
ulation and wind tunnel tests, was done using an unswept, 
symmetric wing model to determine the basic characteristics 
of an active flutter control system. A simple aeroelastic 
model consisting of a cantilever plate with a follower force 
at its free end was analyzed both experimentally and theore- 
tically in attempting to actively control flutter in Refs. 7 
and 8. Analytical analyses of more complicated structural 
models applicable to the B-52 and F-4 respectively are des- 
cribed in Refs. 9 and 10, and in Ref. 11 a control system 
design based on aerodynamic energy is discussed. Each of. 
these analyses, though detailed in many respects, considered 
rather simple sets of feedback parameters in designing their 
control systems. This approach, dictated primarily by sensor 
capabilities, ignores the possibility of using observors or 
estimators in determining unsensed parameters. The analyses 
have, therefore, been severely limited in a number of respects 
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(1) the maximum available increase in flutter velocity with 
a particular set of controllers can be limited, (2) deter- 
mining the gain values that produce this maximum flutter ve- 
locity increase can be tedious for complex systems, (3) se- 
lecting a control system design giving a particular flutter 
velocity increase while optimizing other pertinent design 
parameters is difficult, and (4) design procedures cannot be 
easily adapted to multi-input /multi -output systems. To re- 
move some of these limitations, this thesis incorporates 
three goals: 

(a) Develop and demonstrate organized synthesis techni- 
ques for the design of flutter suppression systems 
using concepts of modern control theory 

(b) Using a control moment gyro as the primary control 
torque generator, indicate and illustrate design 
tradeoffs associated with the general development 
of a practical flutter suppression system 

(c) Establish and identify attractive design features 
of the control moment gyro to warrant its consider- 
ation as a viable alternative to aerodynamic sur- 
faces in producing control forces to suppress flut- 
ter modes . 

1 . 2 Control Moment Gyro 

Although not usually mentioned when possible methods for 
active flutter suppression are considered (4, 12), substantial 
experience has been obtained using control moment gyros (CMGs) 
for momentum exchange in the attitude control of space vehicles 
and satellites (13 - 19) . To familiarize the reader with this 
type of control, the following is a brief description of the 
CMG, including some of the possible advantages to be obtained 
by its incorporation in a flutter suppression system. 

The CMG consists basically of a rotor (or wheel) with 
a known constant magnitude of angular momentum, and a method 
by which the rotor can be precessed. As discussed more thor- 
oughly in Ref. 18, the three basic configurations for the 
CMG are, as shown in Fig. 1-1, the one-degree-of-freedom, the 
two-degrees -of -freedom, and the twin gyro. The one- and two- 
degrees -of -freedom gyros consist of one rotor mounted so that 
it can be precessed about one or two axes respectively. These 
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systems, however, can have undesirable cross -coupling torques. 
The twin gyro consists of two counter-rotating gyro elements 
which are precessed through equal but opposite angles so that 
the control torque is always directed about a single axis. 

In a three-dimensional flutter analysis, the possibility does 
exist that cross-coupling torques occurring for the one degree 
of-freedom gyro could be used to advantage. However, their 
inclusion in the flutter analysis introduces nonlinear terms 
into the equations of motion. Therefore to maintain linearity 
the twin-gyro control system will be used. 

For space vehicle and satellite attitude control appli- 
cations, the control moment gyro has been shown to have ad- 
vantages in power, accuracy, response, and simplicity com- 
pared to reaction -wheel systems (16) . Since output torque 
is proportional to input gimbal rate, the CMG can produce 
large moments on the airfoil at the relatively high flutter 
frequencies without requirements for large angular momenta, 
and hence, large size. The effectiveness is not influenced 
by the complex aerodynamic forces which depend upon Mach 
number and unsteady aerodynamic effects. Consequently, the 
mathematical complexities in the analysis are somewhat reduced 
Furthermore, the inherent lag in the development of forces on 
aerodynamic surfaces does not occur with this system, 
which should result in a reduction in the amount of lead ne- 
cessary for control and higher allowable controller bandwidth. 
Also, deflection angle of the' gyro gimbal is an easily meas- 
ured quantity and can be used directly as a feedback param- 
eter and in estimating other states of the system. 

1 . 3 Thesis Summary 

In Chapter II, the equations of motion for the typical 
section (two -dimensional wing surface) and CMG controller 
are developed. Incompressible flow is assumed and the un- 
steady aerodynamics for arbitrary motion are represented by 
utilizing the Wagner indicial function. The open-loop char- 
acteristics of the system are examined with critical design 
parameters identified. 

In Chapter III, the closed-loop characteristics of the 
system developed in Chapter II are determined. The system is 
placed in standard state variable form with the definition of 
two new states, and the controllability and observability of 
the system are discussed. Full state and reduced state feed- 
back designs for various CMG sizes are determined. 
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b. Two Degrees of Freedom 


/ 



c. Twin Gyro 

FIG. 1-1 CONTROL-MOMENT GYROS 
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In Chapter IV, optimum feedback systems based on a num- 
ber of design criteria are determined. System response to 
gust inputs is also calculated. 

In Chapter V, the equations of motion for a three dimen- 
sional, uniform cantilever wing are developed using assumed 
mode methods. Optimum orientation and location of the CMG 
on the wing are determined. The open loop characteristics 
of the system are calculated and comparisons are made with 
the two dimensional analysis. 

In Chapter VI, the closed loop characteristics of the 
system are examined including a discussion of the effects 
on controllability and observability of increasing the num- 
ber of vibrational modes in the analysis. Critical gains in 
full state feedback systems are identified, and a number of 
reduced state systems are determined. The problem of higher 
mode instabilities induced by the control is analyzed and 
methods for reducing the probability of its occurrence are 
discussed. 

Chapter VII presents a summary of results, conclusions, 
and areas of further research. 
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II. TYPICAL SECTION WITH CMG: 

OPEN LOOP CONTROL CHARACTERISTICS 

2 . 1 Wing-Controller Equations of Motion 

The equations of motion of the typical wing section, 
illustrated in Fig. 2-1, appear in most aeroelastic refer- 
ences including Ref. 20, and will not be derived here. The 
equations of motion for the control -moment gyro are derived 
in Appendix B and their inclusion in the typical section 
equations is relatively straightforward. The resultant equa- 
tions can be written 


m* h + C^h + 

S a a - -L d - L 

(2-1) 

S a h + T-a’a + 

C a a + H w a = Md + M a 

(2-2) 

-HwS + 

1gz ct = T c 

(2-3) 


where the wing structural and geometric parameters, CMG 
parameters, and CMG orientation are defined in Fig. 2.1, the 
list of symbols, and Appendix B, respectively. The spin and 
precession axes of the gyro have been selected so the output 
torque generated by the CMG is about the elastic axis of the 
typical section. The effect of alternative orientations is 
discussed in Chap. V, when the three-dimensional case is con- 
sidered. Structural damping contributions are considered 
negligible » 



FIG. 2-1 TYPICAL SECTION GEOMETRY 


For transient motion and incompressible flow, the ex- 
pressions for aerodynamic lift and moment per unit span are (21) 
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+ 2ffpbU 


f j^h(t) 4- Ua(t) + b(%-a)a(t)J 4>( r - t)dt 

J 0 dt 

M a “ Trpb^j^bah - Ub(^ - a)a - b^ ('/» + a^)aj (2-5) 

+ 2j7pb 2 U(% + a)j" -4r[h(t)+Ua(t)+b(%-a)a(t)J <l>(r-t)dt 


where t = Ut/b is a dimensionless time measured in semi- 
chords traveled after the start of the motion at t = 0. The 
superscript dots, however, denote the derivative with respect 
to physical time, t. (t) is the Wagner indicial admittance 
function describing the circulatory lift buildup after a 
sudden change in incidence, and is discussed in Appendix A. 

Substituting Eqs. 2-4 and 2-5 into 2-1 and 2-2, and then 
dividing Eq. 2-1 by irpb J and Eq. 2-2 by Trpb^ produce the equa- 
tions 

[m + l] £ + H-. jj- + <x a ^ - a)a + £a 


+ 2 


U 

b 


f4| 

X dt 


h + U 



- t)dt = - L D; 

(2-6) 


(x a /A - a) ^ + [r j. + ('/.+ a^)] Q + - a)a 

+ - 2 &<%+«) 

x \k+£ a + (% - *)a] *(r- t)dt + G^a = M D ; (2-7) 

L b b J 

-G 0 w a + o - — . (2-8) 

iGz 

In order to simplify future operations with these equations, 
Eq. 2-6 is multiplied by (% + a) and added to Eq. 2-7, which 
yields 



(2-7a) 


[x a /* + t±(k + a) + h] £ + + a > K 

+ + a) + + %(%-a)]a + j~« + rjo^a + Gi^a = 0. 


Equations 2-6, 2-7a, and 2-8 are then the equations of motion 
for the wing-controller system. 

2 . 2 Open-Loop Characteristics of Wing-Controller System 

Taking the Laplace transform of the equations of motion, 
with the initial conditions and the disturbing force equal 
to zero, gives in matrix form 






f-» — 

S*[ M+l+2 ^ l-2e)$( s) s)j s O 

b 


h(s) 

0 


0 

|+a)+J [x^i( £+a)+r^+$(£~a) }s 2 + ~s+r^u^i G l u a S 


a(s) 

a 

0 

! o <* 2 


°(s) 


V> 

I 

L J 



^ *Cz J 


(2-9) 


where <!>(s) is the Laplace transform of the indicial admit- 
tance function discussed in Appendix A. 

It is of interest to examine the open-loop characteris- 
tics of this system with T c = 0. 

For convenience and to nondimensionalize all the param- 
eters, the substitutions 

_ s 


ba a 

*>( g) = s$(s) = °.5s 2 + 0.281 us + 0.01365 u 2 

s 2 + 0.3455 us + 0.01365 u 2 

2-9, and the characteristic determinant 


are made into Eq. 
becomes 
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(y+l)5 2 +2ucp(s )s+(~) n -a)s 2 +ujl+(l-2 a )tp(g)j; + 2u 2 ?i( 


s) 0 


j*^+M(§+a)+jJ s + ( (j ^)n(i+a) J+a)+r^i+J(i.-a)j s 2 +us+r^i G^s 


-G a 
o 


= o 


(2-10) 


The characteristic equation is written out in Appendix C and 
has the form 


i 2 { A 6 I 6 + u A 5 s5 + (u 2A 41+ A 4 2 + A 4 3 G 2 )s 4 + 5(ii 2 A3i + A 32 + A 3 3 G 2)5 3 


+ [u^A +u 2 A + A +(u 2 A +A )G Is 
L 21 22 23 24 25 2j 

+ u [u 2 A +A +(u 2 A +A )G Is 

L 11 12 13 14 2 J 


( 2 - 11 ) 


+ u^ (u^A +A +A G ) 1 = 0 
01 02 03 2 S 


where the A coefficients are defined in Appendix C and are 
functions of f± t a, Xg, r^, and (wh/o>a) > and G 2 is the gyro 
parameter equal to the product 

0lG o - % . 

n P b ^ l Gz^a 

In order to determine the effect of the CMG on the 
dynamic characteristics of a typical section, including flut- 
ter velocity, the following values were selected for the non- 
dimensional wing parameters (similar to those in Ref. 22) 

M = 10.0 

a = -0.4 

= 0.2 

= 0.25 

0> a y = 0.33 . 
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Substituting these values into Eq. 2-11 yields 


s 2 |24.88s 6 +15.97 us 5 +(0.36u 2 +36.67+11.0 G2)s 4 
+ u(-l, 31 u 2 H-17«84+4.8 G 2 )s 3 
+ [-0.40 u 4 +2.76u 2 +8.25+(0.71u 2 +3.3)G2] s 2 
+ u [-0 . 085u 2 +2 . 85+ (0 . 006u 2 +l , 14 ) ^2 J s 
+ u 2 [-0.009u2+0. 113+0. 045 G 2 ] } 

= 0 . 


( 2 - 12 ) 


Figures 2-2a and 2-2b show a comparison of the root 
loci for this system for three values of G 2 including G2 = 0, 
i.e„, without CMG. The root loci are drawn for variations 
of the velocity parameter u = U/bw a , which will imply, in 
this analysis, variations in U with <u a and b held constant. 
Roots marked "x" indicate u = 0. From the figures, it can 
be seen that for the particular values chosen for the non- 
dimensional parameters, the critical roots associated with 
the onset of flutter correspond to the torsion branch of the 
root locus plot. The frequencies associated with the bending- 
torsion modes approach one another as the flutter velocity is 
approached for each case, which is a result often occurring 
in classical bending-torsion flutter. For increasing values 
of u beyond its flutter value, the frequencies corresponding 
to the bending and torsion branches of the root locus tend 
to converge to the same value. The roots that remain fixed 
at the origin for varying u are associated with presence of 
the CMG in the system. 


For increasing values of G 2 , corresponding increases in 
flutter velocity and flutter frequency occur. This can be 
seen more explicitly in Fig. 2.3 which shows that increasing 
G 2 from 0.0 to 5.0 results in increases for up from 1.65 to 
3.35 and for co p/o> a from 0.80 to 1.11. In general, increases 
in the value of G 2 would correspond to increasing the size of 
the CMG. However, for a constant value of angular momentum, 
increasing G 2 could correspond to a decreasing value of Iq z 
with fixed angular momentum, which would be associated with 
a decreasing CMG size. 
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FIG. 2-3 VARIATION IN VELOCITY AND FREQUENCY AT FLUTTER 
WITH G„ 


In the next chapter, the effect of feedback on this sys- 
tem will be examined. 
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III. FEEDBACK CONTROL OF TWO DIMENSIONAL 
WING- CONTROLLER SYSTEM 


3 . 1 State Variable Formulation 


The characteristic equation developed in the previous 
section is 8th degree, where the two roots located at the 
origin in Figs. 2-2a and 2-2b are associated with the pres- 
ence of the CMG in the system. If feedback of the gimbal 
angle a were not required, the degree of the characteristic 
equation could be reduced to seven. However, passive or 
active feedback of gimbal angle is necessary to maintain the 
validity of the small-gimbal-angle assumption in the equa- 
tions of motion. Thus, the minimum number of states to de- 
fine this system is eight. It would be advantageous for 
sensor and estimator design requirements if it were not nec- 
essary to obtain information on all states to adequately 
suppress the onset of flutter for reasonable increases in 
the velocity. To obtain some insight into which parameters 
are most significant in the control law, an analysis using 
arbitrary dynamics is done in this section. In this method, 
the eigenvalues of the uncontrolled system are shifted to 
arbitrary new locations using state variable feedback. This 
requires that the system be modally controllable with the 
particular set of states chosen to define the system. To 
determine the controllability of the system, the equations 
of motion (Eqs. 2-6, 2-7a, 2-8) are first put into standard 
state-variable matrix form, which can be written in general, 
with no disturbing forces 


x = Fx + G up (3-1) 


where x is the state vector and Up the control vector. With 
^(s)^ N^/D^, the first equation in Eqs. 2-9 is multiplied 

by D^, , and then the inverse Laplace transformation of Eqs. 
2.9 is taken assuming zero initial conditions. The result- 
ing equations can be written 

•« 

B lb" + (Bjdj+Uujjf (B^qU 2 + 2u 2 n 1 + B 4 )|* 



+ u^(B 0 d + d 1 + B^n., + l)d 
z o 1 3 1 


+ u 3 [d 0 + B 3 xi 0 + 2n 1 ]a + 2n Q u 4 a = 0 ; 

[b 2 + Bi(% + a)] |+ B 4 (% + a)£ + [fi 2 (% + a) + 
+ 5a + r^/xa + = 0; 


(3-2a) ' 

4 m - + a2 ] 5 

(3-2b) 


~ G o a + a 


t 2 
a> 

Gz a 


(3-2c) 


Equations 3-2 are equivalent to Eqs . 2-6, 2 -7a, and 2-8 in- 
sofar as they have the same characteristic equation. For 
certain nonzero initial conditions, the dynamic response of 
the two systems would be different. 

Before writing Eqs. 3-2 in standard form, a decision . 
must be made as to which quantities are considered states. 
There are many possible realizations of the system of equa- 
tions, and the procedure here will be to select those states 
which would be the least complicated to individually sense, 
and can be explicitly identified. The following quantities 
are therefore selected as the states of the system: 


hi = h 

(3-3) 

h 2 “ 

(3-4) 

. a 

(3-5) 

“2 * “1 

(3-6) 

CT^ = <7 

(3-7) 

ct 2 " *1 . 

(3-8) 

- (B i ir + b 2 V 

(3-9) 

* h 

^1 + (B^d^ l)u 2 l . + u(B 2 di + 1 + 6^2)02 

(3-10) 
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Substituting Eqs. 3-3 to 3-10 into Eqs. 3-2 and writing in 
matrix form yields Eq. 3-11, where the F-jj are defined in 
Appendix Cs 



A necessary and sufficient condition for arbitrary re- 
location of the system eigenvalues by state-variable feed- 
back (i.e., modal controllability) is that the so-called con- 
trollability matrix be nonsingular. The controllability ma- 
trix is defined to be, 

C — I" G | FG ! - - • ! F n_1 G] . (3-12) 
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where n is the number of states defining the system. It can 
be shown that this matrix is nonsingular for the parameters 
selected in this analysis with G^, G Q >0, u > 0. 

Thus, using only one CMG -controller and full state feed- 
back, the flutter velocity can theoretically be increased to 
any value by just making all the real parts of the system 
eigenvalues have negative values. This generally, is not 
only true for systems utilizing a CMG controller, but also 
can be shown to be the case if a single aerodynamic surface 
were to be used. This is contrary to implications given in 
Ref. 11, in which Nissim states that the typical section 
system "has two degrees of freedom, translation (bending) 
and rotation (torsion). Therefore, two control surfaces are 
required in order to maintain precise control of the two de- 
grees of freedom of the main surface." However, implicit in 
that requirement, although it is not entirely clear what is 
meant by precise, is that the number of feedback parameters 
is limited to only four quantities, a, a, h, and h. 

3 . 3 Full-State Feedback Gain Determination 

Another method for determining if the system is modally 
controllable is to actually solve for the feedback gains in 
the control law, which relocate the eigenvalues to selected 
locations. Obviously, if the gains can be determined, the 
system is modally controllable. Two methods will be given 
for the determination of the required feedback gains. The 
first is rather straightforward and can be generally more 
convenient to use on low order systems, since the evalua- 
tion of the controllability matrix is not required, and ac- 
tually the system equations need not be arranged in state 
variable form. The second method is mathematically more in- 
volved, but provides a generalized format for more conve- 
niently calculating the gains in a large order system. The 
description of this method is based on the procedure out- • 
lined in Ref. 25. The control law will be assumed to be of 
the general form 




or in vector notation 


T 

Uc = -K x 


I Gz <u a 


- - [-K x ,^ 2 , -K 3 , -K 4 , +K 5 ,+K 6 , -K 7 , -K g J 


h l 

h 2 

«1 

a 2 

°l 

a 2 

€i 

\J 2 j 


(3-14) 


The basic problem of selecting gains to arbitrarily 
place the system eigenvalues can be stated as follows: 

Given the single-input system r 


x = Fx + gu c 


(3-15) 


(where for application to Eq. 3-11, Uc is a scalar and £ is 
a vector) with the characteristic polynomial for uq = 0, 

X(s) - det (sI-F) = s n +f x s 11 " 1 + f n , (3-16) 

find the gain vector K such that the characteristic polyno- 
mial of the system 

x = (F - gK T )x (3-17) 

has the specified form 

*k( s ) = s n + g 1 s n_1 g n . (3-18) 

Whenever the system contains a small number of states, 
a solution for the gains can sometimes be generated more 
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conveniently by directly expanding the characteristic deter- 
minant of Eq. 3-16 and equating the coefficients of like 
powers of s with those of the desired characteristic equa- 
tion, Eq. 3-18. For example, for the particular system des- 
cribed in Eq. 3-11, a solution can be developed by substituting 
the Laplace transform of Eq. 3-13 into the transform of Eq. 

3-2 which results in the following characteristic determi- 
nant : 


B^ 4 + uCb^-H)! 3 + [(B l d o +2n 1 )u 2 +B <l ]s 2 B^ 4 + ( B^+UB^ ^ s3 


+ ( 2 u 3 n +B d,u)s + B.d Q 
o 4 X 4 o 


. 2 


[B^ + Bjfj + a)l s 2 l + a) 

-K 8 B iS 3 + £-X 7 B 1 -K fl fi(B 1 <J i+ l)ji 2 


2.2 


+ (B-d +d.+B,.n +l)u s 
v 2 o 1 3 1 

3 4 

+ (d +B n +2n t ) a S + 2n 5 
x o 3 o 1 o 

[B 2 (J+a)+r 2 p-i^a 2 ]5 2 + us + r^i GjS 
-K g B 2 3® + [-K 7 B 2 -K 8 0(8^ 

♦ 1+B 3 n 2 )]s 2 I 2 -« 4 s* 5 

■- Cv°.)» " K 3 


= X. 


K 


(3-19) 


Expanding this determinant and equating coefficients estab- 
lishes eight equations, which can be solved simultaneously 
for the eight gains to arbitrarily relocate the eigenvalues 
of the system. 


For large order systems the above approach could be 
tedious, and a procedure more amenable to computerized solu- 
tion is needed. The following method is presented from Ref. 
25. From Eq. 3-17 

*K(s) = det (si - F + & K T ) 

- det (si - F) (I + (si - F)' 1 & K t ) 

- det (si - F) 'det (I + (si - F^&K 1 ) 

r X(s)‘ (1 + K t (sI - F)" 1 g.) (3-20) 
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Hence , 

Y K (s) - X (s) - X(s) • K T (sI - F)- 1 & 

= Yi 9 11 " 1 + Fn-1 S + Yn (3-21) 

where 


-1 y l s ”" 1 + y n-l s + F n 

x rr — ~ " 


K 1 (si - F r x S. = 


s n + f l S 

Now, as can be directly verified. 


TTT 


L n 


n-3 


(si - F)" 1 = 1 [s’ 1-1 + (F + fll)s n ’ 2 + (F 2 + fiF + f 2 l)s 

X(s) L 


+ (jP -1 + f + • * • fn -il)j 


(3-22) 

Substituting Eq. 3-22 into Eq. 3-21 and equating coefficients 
of like powers of s gives the following system of equations: 

*1 - k t & 

>2 = K T Fg. + f]K T g. 

Y 3 = K T F 2 g. + fiK T Fg. + f2K T a 


Vn = -I- f iK T F n ' 2 g. + £n-lK T g. (3-23) 

This can be written in matrix form as 
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- rc t k 


K T £. 


s. T 

k t f& 


R T F T 


= R 


K T F n_ lg_ 


gT p T (n-1) 


(3-24) 


where C is the controllability matrix, and _y is the vector 
whose elements are the differences in the coefficients of 
the characteristic polynomials defined in equations 3-16 
and 3-18, i.e.. 


g l - f l 
§2 " f 2 


g. 


n 



(3-25) 


Thus, the gain vector can now be obtained from Eq. 3-24: 


K * C'V 1 ! (3-26) 

This will have a solution, since R is clearly non-singu- 
lar and, for the controllable system that is considered here, 
C is non -singular. 
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з . 4 Design by Arbitrary Dynamics (Pole Placement) 

An extensive parametric study was done, utilizing equa- 
tions obtained from Eq. 3-19, in order to determine some of the 
characteristics of the wing-controller system. Variations in 

и, G2, and root locations were considered. Although later de- 
sign specifications may dictate additional requirements on the 
root locations, the one stability criterion established for 
arbitrary root placement was that the minimum negative value of 
the real part of any root would be greater than 0.05<u a when u 
is greater than its critical value for the wing without the 
CMG (i„ e. , u >1.65). Consequently, two particular situations 
arose in the analysis, depending on the values of G2 and u: 

(1) one conjugate pair of roots had positive real parts or did 
not meet the stability criterion, and thus, to maintain an 
adequate stability level, these roots plus the two roots at 
the origin (Fig. 2-2) had to be relocated by proper selection 
of the control gains; (2) all roots except those at the origin 
met the stability criterion, with no positive real roots. The 
roots that are relocated in order to meet the stability crit- 
erion will be referred to as the critical roots. When judging 
the merits of various root locations for particular values of 
G2 and u, the relative size of each set of control gains was 
the basis of comparison, with smaller gains implying a more 
desireable design. Though not a direct indication of the qual- 
ity of a particular design, the size of the gains do provide 
a measure of the relative minimum values of the torque and 
power requirements among designs constrained to the same stab- 
ility criterion. This is verified by the results of Chap. 4. 
Low gain values also facilitate the determination of adequate 
reduced feedback systems, since elimination of certain feed- 
backs as unnecessary can be made when their values are negli- 
gible. The following results and conclusions were obtained 
from the parametric analysis 

a. The lowest gain values are obtained when only the crit- 
ical roots are relocated and are placed as near to the 
imaginary axis as allowed by the stability criterion. 

b. When relocating only the roots at the origin, the 
values of the gains are reduced by placing the roots 
on the real axis. 

c. When there are four critical roots, the gain values 
can be reduced (1) by proper adjustment of the imag- 
inary part of roots originally located at the origin, 
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and (2) by relocating the torsion-branch roots at 
approximately the same distance from the origin as 
they were with no control- 

d. Gain scheduling, in which the gains are made a 
function of freestream velocity, can be used to 
reduce the gain values over a range of veloc- 
ities. Depending on the value of G2> all the gains 
need not be scheduled, but even with reduced sched- 
uling, the added complexity may be prohibitive for 
a particular design. The following three examples 
illustrate the effects of gain scheduling. Since 
the wing is stable without the CMG for u <1.65, 
only the feedback characteristics for u >1.65 were 
considered of interest. 

Table 3-1 shows, as a function of velocity and for 
three values of G 2 , the gains needed to place the 
eigenvalues of the system at the desired new loca- 
tions as shown in Figs. 3-1, 3-2, and 3-3. The 
root locations were chosen to produce minimum gain 
values using the guidelines indicated above in a, 
b, and c. In the characteristic equation, all the 
gains except K 5 and Kg appear multiplied by the 
gyro term, G^ . Thus in Table 3-1, = G^K^ is 

used as a convenient parameter to describe the 
feedback characteristics, since, for a given G 2 
and set of root locations, will remain constant 
for variations in G^ and G 0 . Figures 3-1 to 3-3 
show three different control situations: 

(1) G 2 = 0.221: There are four critical roots 

for all u. 

(2) G 2 = 3.8: For each u except u = 3.0, only 

two critical roots exist. 

(3) G 2 = 2.12: For u<2.4 there are two critical 

roots, and for u > 2.4 there are four critical 
roots . 

For G 2 = 3.8, the gains remain relatively constant 
over the range of velocities considered except at 
u = 3.0. This indicates that a set of constant 
gains could be selected which would place the roots 
in approximately the same locations. For G 2 = 0.221, 
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Table 3-1 


SCHEDULED GAINS FOR FLUTTER SUPPRESSION TO u = 3.0 


u 

Kl 

K 2 

Ki = G x 
K 3 K 4 

K i 

*5 


k 7 

Kg 


— 


( 1 ) g 2 - 

0.221 




1.75 

.043 

.070 

-.110 

-.060 

.148 

.207 

-.021 

-.002 

2.0 

.062 

.118 

-.203 

-.105 

.231 

.340 

-.030 

-.016 

2.3 

.082 

.166 

-.195 

-.000 

.291 

.456 

-.036 

-.018 

2.5 

.098 

.198 

-.140 

- 0.0 

.317 

.522 

-.039 

-.020 

2.7 

.115 

.236 

-.086 

.05 

.332 

.582 

-.061 

-.021 

3.0 

.127 

.317 

-.070 

.06 

.300 

.665 

-.042 

-.022 




( 2 ) G 2 = 

2.12 




1.9 

.014 

.025 

.231 

.032 

.005 

.10 

0.0 

.006 

2.1 

.014 

.024 

.235 

.035 

.006 

.10 

0.0 

.005 

2.3 

.015 

.023 

.241 

.040 

.006 

.10 

0.0 

.005 

2.5 

.015 

.020 

.035 

-.265 

.024 

.2175 

-.0164 

-.0201 

2.7 

.091 

.163 

.056 

-.300 

.071 

.339 

-.024 

- .022 

3.0 

.145 

.315 

.109 

-.360 

.105 

.490 

-.026 

-.026 




( 3 ) G 2 = 

3.80 




1.9 

.026 

.045 

.415 

.058 

.008 

.10 

0.0 

.011 

2.1 

.026 

.042 

.422 

.063 

.008 

.10 

0.0 

.010 

2.3 

.027 

.041 

.431 

.071 

.009 

.10 

0.0 

.009 

2.5 

.029 

.042 

.445 

.082 

.010 

.10 

0.0 

.009 

2.7 

.033 

.046 

.464 

.098 

.012 

.10 

0.0 

.009 

3.0 

.010 

.013 

.257 

-.285 

.016 

.197 

-.011 

-.0155 
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Imaginary 





FIG. 3-3 ROOT LOCATIONS FOR G_ = 3.8 



each of the gains varies somewhat, with only feed- 
backs of h and a having a significant variation 
with velocity. For G 2 = 2.12, each gain does show 
significant variation with velocity, caused pri- 
marily by the change from two to four critical 
roots as the value of u increases beyond 2.4. 

Figures 3-4, 3-5, and 3-6 show a comparison of the 
root loci for G 2 = 0.221, 2.12, and 3.8 respectively, 
using constant gains and gain scheduling. Figure 
3-6 shows that constant gains can be used with very 
little change in the root locus. Figure 3-4, G2 = 
0.221, shows that there is very little change in 
the noncritical roots using constant gains. The 
critical roots become located on either side of 
their scheduled value. For G 2 = 2.12, Fig. 3-5 
shows what differences may result when using con- 
stant gains. The roots corresponding to the torsion 
branch of the root locus with constant gains follow 
the bending branch associated with scheduled gains, 
and the bending branch of the root locus with con- 
stant gains results in the critical flutter condi- 
tion. 

Thus, it appears that to obtain the best results 
using constant gains, the value of G 2 should be 
selected to give an open-loop flutter velocity 
greater than the maximum flutter velocity desired 
for a particular application. 

e. Relaxing the stability criterion, such that the 
critical roots can be placed nearer the imaginary 
axis, can result in reduced values for the gains. 

The size of the gains, however, may not provide an 
accurate measure of the relative values of the torque 
and power among designs constrained to different 
stability criterion. 


f. Although feedback of eight states is necessary to 
place exactly the set of roots at a desired loca- 
tion for particular values of velocity and G 2 , 
adequate stability can be achieved with reduced- 
feedback systems. A number of possible designs 
using various combinations of the feedback param- 
eters can be developed, depending on the values of 
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u, G 2 and the eigenvalues of the system. One ex- 
ample of a reduced-feedback system can be obtained 
from Table 3-1 for G 2 = 3,8. The significant gains 
correspond to feedback of a and a , and a design 
using only these two feedbacks, along with o to null 
the gyro, results in a root locus very similar to 
the eigenvalues shown in Fig. 3-6. 

The simplest feedback system for any value of G 2 
can be determined by examining the characteristic 
determinant, Eq. 3-19. It can be seen that feed- 
back of a(K 4 ) has the same influence on the roots 
in the characteristic equation as G 0 . Consequently, 
increasing the value of the gain parameter K 4 is 
equivalent to increasing G 2 . Figure 3-7 shows 
the root locus for a system using only two feed- 
backs characterized by the gains 


(K 4 + G 2 ) - 3.8 

k 5 = 0.10 

For G 2 = 3.8, this system reduces to one requiring 
the feedback of only a . 

A comparison of the response characteristics of this 
system with systems using all eight states as feed- 
backs will be done in Chap. 4. 

3.5 System Observability 

In the preceding analysis, total knowledge of all the 
necessary states for use in the control laws was assumed. 

This information could be obtained through sensor measure- 
ments or through a properly designed estimator if it were 
not possible or convenient to measure directly the desired 
states. The sensor problem would be reduced to a minimum if 
all the states could be estimated from one measurement, the 
most convenient of which would be o . If the state of a sys- 
tem at any time can be uniquely determined by the available 
measurements, the system is considered observable. The ob- 
servability of the wing-controller system, assuming only a 
is measured, can be qualitatively determined by examining 
Eqs. 3-2. Mathematically, observability can be determined 
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by the requirement that the observability matrix be non- 
singular. The observability matrix is defined to be 


r 



HF 


0 = 


HF 


n-1 


(3-27) 


where H is the output matrix defined from 


2. = Hx 


(3-28) 


and y. is the measured quantity. With only a measured, 

H = [o 000100 Oj . (3-29) 

It can be shown that the observability matrix is non- 
singular for the parameters selected in this analysis for 
G 0 , Gi >0, u > 0. 

Thus, it appears a may be the only quantity necessary 
to measure, since stability can be achieved with only o feed- 
back. If more control of the eigenvalue locations is desired, 
the other states can be estimated from this measurement. 
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IV. SYSTEM DESIGN USING OPTIMAL 
CONTROL 


4.1 Optimal Control Using Quadratic Synthesis 

In Section 3.4, the relative values of the gains in 
the control law were used as a basis for comparing various 
eigenvalue locations. A criterion more directly associated 
with system performance, which results in stable, linear 
feedback control, can be obtained by minimizing a cost 
function made up of integral quadratic forms in the state 
and control. Determining the control requirements in this 
case is basically similar to the optimal regulator problem 
as discussed in Ref. 26. The following is a brief descrip- 
tion of this method. 

As has been shown previously, the equations of motion 
can be written. in the form 

x = Fx + Gu c> (4-1) 

It is desired to minimize a cost function of the form 


J = - ( (x^Ax + u c ^Bu c )dt (4-2) 

2 

where superscript T indicates the matrix transpose, and A 
and B are symmetric weighting matrices selected to obtain 
"acceptable" levels of x and u,, . The control u c which 
minimizes J can be obtained by solving Eq. 4-1 simultan- 
eously with the Euler -Lagrange equations 


11 

H , 
• *< 1 

3x 

(4-3) 

aH = 
dUc 

0 

(4-4) 


where 

y4 - %x^Ax + ^Ur^Bu c + A T (Fx + Gu*-,) 
and A is a vector of Lagrange multipliers . 
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(4=5) 


Substituting for ^ in Eqs, 4-3 and 4-4 yields 

A - - Ax - F T A 

u c = - B _1 G T A.. (4-6) 

Substituting Eq. 4-6 into 4-1 and combining with Eq. 4-5 gives 


— m 

X 


’ F 

-GB -1 G T ' 


x" 

■ 


-A 

-ft 


A 


One method of solving Eq. 4-7 is by using the substitution 
\ - Sx, where S is a matrix. This results for the regula- 
tor problem in determining the steady-state solution of a 
matrix Riccati equation for S of the form, 

SF + F T S - S T GB-lG T S + A = 0 (4-8) 

The feedback gain matrix, K, in the control law u c = -Kx 
is then obtained from 

K = B-!g t S. (4-9) 

Another method of solution is by eigenvalue decomposition as 
described in Ref. 27. 

Obviously, critical to any solution of this problem is 
the choice of weighting matrices. One "rule of thumb" in 
selecting A and B is to choose diagonal matrices whose 
elements are equal to 

1 and 1 

( x imax)^ > (ucmax)^ 

respectively, where xi max and u Cmax are the maximum 
"acceptable" values of xj[, the elements of state vector, 
and U£ respectively. It must be remembered, of course, 
that the weighting matrices are only relative, and the actual 
maximum values and their relative sizes occurring in practice 
will depend on the system response to initial conditions and 
disturbances . 
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As a first step in this introductory analysis, the 
diagonal elements of A are considered to be very small 
with respect to B, so that no restrictions are placed on 
the state behavior other than stability, and therefore the 
design that essentially minimizes J ' lie dt i s being 
determined. The matrix A can be positive semi-definite 
to the extent that it cannot be identically zero, since, as 
can be seen from Eq. 4-2, the solution degenerates to the 
case, u c = 0 . 


Table 4-1 shows the feedback gains for G2 = 0.221 as a 
function of velocity, resulting from a solution of the 
optimal control problem. 


A root locus of the critical roots is shown in Fig. 4-1. 
The noncritical roots, as might have been expected, are left 
unchanged in the controlled case, and are not shown in the 
figure. 


TABLE 4-1 


GAINS FROM OPTIMAL CONTROL 
ANALYSIS G? = 0.221 
Ki = G 1 K ± 


u 

K 1 

k 2 

k 3 

k 4 

K5 

K6 

*7 

k 8 

1.75 

-.005 

-.010 

-.028 

-.039 

.001 

.057 

-.003 

-.006 

2.3 

-.166 

-.026 

-1.26 

-1.22 

.001 

.555 

-.038 

-.168 

2.5 

-.202 

.052 

-1.56 

-1.45 

.001 

.687 

-.037 

-.182 

2.7 

-.232 

.136 

-1.81 

-1.63 

.001 

.807 

-.037 

-.187 

3.0 

-.272 

.260 

-2.13 

-1.84 

.001 

.971 

-.039 

-.186 


For a particular value of u, it can be seen that if 
the system instability is due to eigenvalues equal to R + jl, 
then to minimize J, the roots in the controlled case should 
be placed at -R + jl. Thus, as the value of u increases, 
the real root of the critical eigenvalue should be placed 
increasingly farther into the left-half plane. Similar re- 
sults occur also for larger values of G 2 and can be re- 
peated on a simpler level by looking at a simple second 
order system with negative damping. For example, consider 


- 33 - 





FIG. 4-1 ROOT LOCUS FROM OPTIMAL CONTROL ANALYSIS 
C 2 = 0.221 
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x - 2Rx +c j^x = Uc . 
For Ug = 0, the eigenvalues are 

s = R + jl 


with 

I - 

To minimize 

J = I (x-^Ax + bu£)dt 
with the diagonal elements of A small, 

u r = -4Rx. 

The resulting eigenvalues are 

s - -R + jl. 

Since six of the eight foots remain essentially unchanged 
in the higher -order system when going from the uncontrolled 
to the controlled case, this system takes on characteristics 
similar to those of a simple harmonic system. 

The large gain values appearing in Table 4-1 as com- 
pared to those in Table 3-1 are in part due to the tighter 
control requirements resulting from moving the pair of 
critical roots farther into the left-half plane. Primarily, 
however, the large values are due to the relocation of the 
roots originally at the origin. It was found in Section 3.4 
that proper placement of these roots can result in lower 
gains (and for the particular value of G 2 used in this 
analysis, large gains result for placement of these roots 
near the real axis). Thus, whether gain values can be used 
to indicate a desirable design as was done in Section 3.4 
depends on what particular function is being optimized. 

Because of the placement of the critical roots in 
achieving the minimum integrated value of control torque, 
the initial values of torque tend to be large. This feature 
of the design can be partially alleviated by redesigning the 
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control for larger values of 02° In fact, for practical 
design considerations, the maximum values of torque and 
power may be more critical parameters than the time inte- 
gral of control torque. The effects on these parameters 
and other response characteristics resulting from variations 
in u, G 2 » and root locations are investigated in Section 
4.2. 


4.2 System Response to Gust Impulse 


The aerodynamic lift and moment produced by gust-like 
disturbances in the airstream are a major source of external 
loads on a lifting surface. The lift and moment distur- 
bances created by an impulsive upward gust striking the 
leading edge at time t = 0 can be described by (20) 


Ld * 27rpbu2 


r 


d 

dt 


Wg . 

u (t) 


Mr 


= 2npb^(k + a)U 2 


f i. 

\m] 

K dt 

_U 


* ( »• - t)dt (4-10) 

V' ( T - t )dt (4-11) 


where ^g = WgS(t) is the impulsive gust velocity, and ^(r) 
is the Kussner indicial function for a sharp-edged gust. 
Substituting for Lp and Mp in Eqs . 2-6 and 2-7 results in 


i 1 

T= 

+ 

x ] l + % ^ 

+ (x af i 


+ 2 2 f I 

b ] dt 

J o 1 

’ h + U 
. h b 

= 

- 2 2* CjL.\ \ 
b 2 ) 0 dt L u 

(t) i 

( x a^ 

- a) Jy + r^fjL + Vb + a : 


“ 2 if (% + a) 

ftp 


U 

b 


i + (% - a) 


aj<t>( r 


-t)dt 


(4-12) 


® + jf ~ a )<* + ^<°a r a 2 ° 




r - t)dt 
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2 4 

b 2 

ft + a) 

1 d£ b 

(t)J ^( r - t)dt ; 

(4- 

-13) 


■ G o“ a i 

+ a = 

Tc . 
!Gz 

(4- 

■14) 


As before, adding the product of Eq. 4-13 and (^+a) times 
Eq. 4-12, rearranging terms, and taking the Laplace trans- 
form gives 

2 , 


- 

(M 


A + l)s 2 + 2 u<p(s)s + H - 


o> 


a 


h(s) 

b 


t 


+ | (x^/i - a)s 2 + |b + (l-2a)<p(s)Jus + 2 ^(s)u 2 |a(s) 

-2 Wg-(s) 

= -2u — *(s); ' (4-15) 

x a^ + + a ) + %Js 2 + (^r'j 4 (% + 

+ | [ x q 4 (% fa ) + r aM + %(%" a ) J s 2 + us + r 2 ^ja( . 
+ Gisa(s) = 0; (4-16) 


-G sa(s) + s 2 a(s) = 


T c (5) 

!Gz <^ a 2 


(4-17) 


Multiplying Eq. 4-15 by and substituting the general 
control law, Eq. 3-13 gives in matrix form 
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n - ^ , r.( n a j.i x ([) (j u^-u2n u^+D )s 2 B s 4 + (B d +1+B n )us 

™v ^l“l ' *' 5 ' V 1 o ' 1 4' 2 2.1 3 2 


v --3 


+ (2u 3 n o + B^d^u) s + B^iT 


[B 2 + B 1 (J + a)]s + B 4 (J + a) 
‘K 6 B 1 5 3 + [-K. 7 D 1 -K ft Q(B 1 d l+ l)] 5 2 

-V - K 1 


h( s) 


b 

a(s) 

a( s) 


- 2 *V s) 


o 

0 


+ (B 2 d 0 +^ 1 +B 3 n 1 +l ) ii 2 S 2 

U . _3. .4 

+B^n +2n, )u s+2n u 
o 3 o 1 o 

2 2 ° 2 
[B 2 (J+a)+r^i .+ £+a js'+us+r^i 

-K 3 B 2 s 3 + [-K 7 B 2 -K 8 .S(B 2 d 1+ l 5 2 + K 4 s + K 5 

+ B 3 n 2 )]s 2 -(K 4 +G o )i-K 3 


(4-18) 


The Laplace transforms of the solutions for a(t), h(t), and 
a(t) are 


a(s) = ^!Is^^(s)D v ^B 2 + Bi (% + a) + KgBijs 
+{k7Bi + Kg (Bid 1 + l)u 

+ k 6^ b 2 + B l(^ + a )]}^ 

+{b 4 (% + a) + K 5 ^B 2 + B X (% + a)j + K 2 } 

+ ^K 6 B 4 (% + a) + KiJf + K 5 B 4 (% + a)^ ; 


^ s " 2 


(4-19) 


-38 



o(s) = - 2u 2 / W g^\ d «^(s)^-K 8 A 6 s 5 - r K 7 A 6 + ( A 5D W l A 6^2 A 5N^8l s 4 

*TgGl\ « j V 1 J 

^2+Bi (lr+-a)J (K 4 +G2) + K 8 [b 2 B 4 (%fa) -u 2 (B^pfl) -r 2 M B x j 


+ 
-K2 
+ 


B2(^a)+rJ^+a 2 j -^7B X Js 3 

j^B 2 +B 1 (%4-a)J iC 3 + |B 2 B 4 (%+a)-r 2 ^B 1 jK- 
- K 2 u-Ki^B 2 (%+a)+r^+i+a 2 j 

uK 8 ^B 4 (V'-a)(B 2 d 1 +l+B 3 n2)-r 2 /A(B 1 d x +l)j J 


-2 

s 


+ |B 4 (V»-a)(K 4 +G2)--uKl-r^K2 js 
+ | B4(^+a)K3~r^/i-K^ | \ ; 


(4-20) 


h(s) - -2u i V 5) D„*<i) ( j rB 2 (%+a)+r^+5+a 2 ] +K 8 b| 

b X& U \ ' L 

+ j K 6 ^ B 2 (%+a)+r 2 ^ +g+a 2 J*+ 

+ ^K 7 B 2 +K 8 u(B 2 d x +l+B 3 n 2 ) 

+ j r 2 ^+K 6 i^K 5 j^B 2 (%+a)r 2 ^+a 2 j +(K 4 +G 2 ) 


: s - 4 


s 3 


+ jK 6 r2 AH-K 5 u +K 3 |s+K 5 r 2 ^; 


(4-21) 


where Y K denotes the characteristic equation. like<*>(t), 

is a transcendental function and must be approximated so that 
*(s) can be determined and Eqs . 4-19, 4-20, and 4-21 can be 
inverted. The most familiar form is 


*(t) - 1 - 0.5e“* 13t - 0.50« 


-t 


(4-22) 


Its Laplace transform is 

f(s) = 0.565us + 0. 13u 2 (4-23) 

(s+0 . 13u) (s+u) 

The torque required for control, T c is determined from 
inverting Eq. 4-17, and the power used by the gyro torquer is 
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simply the product of control torque times gimbal rate, b. 
Figures 4-2 and 4-3 show the control torque and power res- 
ponse to an impulsive gust striking the leading edge._ Two 
convenient non-dimensional parameters are presented, T c and 
p (see List of Symbols). For three values of G 2 , 0.221, 
2.12, and 3.8, the feedback gains used in this analysis 
correspond to the values in Table 3-1 for u = 2.5, and the 
value of G-^ is assumed constant at 3.5. 

As can be seen from the figures, the maximum values of 
the torque and power parameters tend to decrease with in- 
creasing G 2 . The control torque and power are damped much 
faster for Go = 3,8, where the real root associated with 
the torsion branch of the root locus is farther in the 
left-half plane (Fig, 3-3). 


The response of the gyro gimbal angle to a gust impulse 
is illustrated in Fig. 4-4. For the same set of roots, the 
maximum values of a for each value of G 9 are approximately 
the same. These maximum values, which equal approximately 
0.2 for the three cases considered, indicate the possibility 
(for large -values of w £ /b corresponding to large gusts or 
small wing chords) that the small-angle approximation intro- 
duced in Appendix B to linearize the equations may become 
invalid. For example, G-^ = 3.5, G 2 = 2.12, w / b = 3.0 re- 
sults in 


<j 

max 


33 deg. 


To obtain the effect of varying G 2 on maximum power re- 
quirements, a parametric optimization study was done in which 
the root locations corresponding to the smallest absolute 
value of maximum power for varying values of G 2 and u were 
calculated. For each G 2 and u the minimum power values occur- 
red for approximately those root locations associated with 
minimum gain values. The results are shown in Fig. 4-5. 

Small values of G 2 correspond to situation (1) as discussed in 
Section 3.4 in which four critical roots exist. The large 
values of control torque are needed not only to null the gyro 
gimbal, but to provide energy to increase the flutter velocity 
of the system. For large values of G 2 , there are only two 
critical roots associated with the gyro gimbal angle. The min- 
imum value of the power parameter, P, for a particular u occurs 
,at the transition between two and four critical roots. This 
minimum corresponds to the actual minimum power required if 
lGz<% remains constant as G 2 is varied. 


To obtain an actual estimate of power required, the 
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FIG. 4-4 GIMBAL ANGLE RESPONSE TO GUST IMPULSE 


following set of values was selected, compatible with the 
chosen values of and G 2 : 

TABLE 4-2 


G 2 

= 3.8 

H 

= 6.0 k s - m2 

w 

sec 

' 

l Gz 

= 0.09 kg-m^ 

"a 

= 10 Hz 

b 

= 1.5 ft. 


The corresponding maximum value of power is 0.007 HP. 

Also shown in Fig. 4-5 are the values of power associated 
with two of the reduced feedback cases discussed in Sec. 3-4, 
and identified on the figure. These values indicate that by pro- 
per selection of the root locations a reduction in feedback 
complexity can be made without a significant increase in power 
requirements (illustrated byO). However, since the system ei- 
genvalues cannot be arbitrarily determined with reduced feed- 
back systems, the noncritical root locations and the stability 
level of the controlled system may be sufficiently altered 
to produce increased power requirements. The maximum gimbal 
angles corresponding to the root locations used to determine 
Fig. 4-5 are shown in Fig. 4-6. The minimum values of a for 
each u occur at the value of G 2 corresponding to the minimum 
value of P in Fig. 4-5. These two curves indicate that selecting 
a value of G 2 corresponding to a value of u^ equal to the maxi- 
mum desired flutter velocity results in the "best" design, with 
regard to minimizing values of power, while keeping variations 
of a to a minimum. This, of course, would be subject to size 
constraints on the CMG itself. The large values of <r/(w g /b) 
may be prohibitive for any practical design. These values can 
be reduced by relocating the critical and noncritical roots of 
the system with corresponding increases in power requirements. 
From Eq. 4-20, it can be seen that increases in Gi, holding G 2 
constant, can also result in a decrease in cr. These would be 
some of the many tradeoffs to be made in the more detailed design 
procedure used when developing an actual system. 
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FIG. 4-6 MAXIMUM GIMRAL ANGLE CORRESPONDING TO THE MINIMIZED 
MAXIMUM POWER RESPONSE 
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v, OPEN -LOOP CONTROL CHAR lACTERISTICS 
OF A THREE-DIMENSIONAL CANTILEVER WING WITH 
CMG CONTROLLER 


5 . I CMG Orientation 


The basic approach to the design analysis of the three 
dimensional wing parallels the typical section analysis 
developed in Chap. II. The wing geometry is shown in 
Fig. 5-1. It is assumed to be an unswept, constant chord, 
cantilever wing, with rigid chordwise sections, having a 
straight elastic axis and uniform structural properties. 
Before developing the equations of motion for the wing con- 
troller system, the most effective orientation of the CMG 
should be determined. As can be seen from Figure 5-2, six 
possible orientations exist, depending on the direction of 
the gyro's precession and spin axes (for illustrative pur- 
poses, only a single rotor of the twin-gyro controller is 
shown). To be most effective, the output torque produced 
by the CMG on the wing, Mq/ w , should be directed about an 
axis of rotation of the wing (i.e., the X or Y axis). 

Thus, cases 5 -2b and 5-2e, which provide moments about the 
Z-axis, would obviously not be a suitable choice. A pre- 
liminary analysis has also shown, for the particular set 
of parameters selected in this thesis (see Section 2.2), 
that directing the CMG output torque about the X-axis, 
cases 5-2c and 5-2f, does not provide adequate control 
capability. Consequently, cases 5-2a and 5-2d appear to be 
the most suitable CMG orientations. For low control torque 
levels, the difference between the two orientations can be 
considered negligible, and for the following analysis case 
5 -2a was selected. The equations of motion for the control 
moment gyro in this orientation are derived in Appendix B. 

5 . 2 Wing-Controller Equations of Motion 

The dynamical equations of motion for the coupled tor- 
sion and bending of a cantilever wing are derived in a number 
of classic aeroelastic references, including Refs. 20 and 21. 
These equations, with the inclusion of the CMG terms from 
Appendix B, can be written 


mh(y, t)-hnch(y, t) S(y-y G )+EI^|^-^d-S a a(y,t) 


+s aG3(y> t)s(y-yG) 
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= -l d -L 


(5-1) 






(5-2) 


S a h(y,t)+S aG h(y»t)s(y-y G )+I a a(y,t)+I aG a(y,t)S(y-yG) 

- + H w <j(y, t) S(y-y G ) = Mq + M a 

- Hwa(y G ,t)+I Gz a(y G> t) « T c (y G ,t) 

With the boundary conditions of a cantilever beam: 

at y = o, h = §£. = a = 0 

3 * ay 

at y =j2 , = ®SL = 0 

3y2 3y3 3y 

The CMG is treated as a concentrated mass and is assumed 
to generate a concentrated torque about the Y-axis. Rota- 
tional inertia terms about the X-axis and CMG effects on 
flexural and torsional rigidity are considered negligible. 

A number of methods are available for developing 
solutions to the above equations. Runyan and Watkins (28) 
developed a solution for a uniform wing with arbitrarily 
placed, concentrated masses by extending the treatment of 
Goland (29), in which the differential equations are attacked 
directly. To obtain flutter information, this is the most 
accurate approach for a uniform wing, but it does not allow 
for a suitable means to determine control gains or to deal 
with wings having non-uniform properties. A practical method 
for performing flutter calculations is based on the assump- 
tion that the motion of the system can be approximated by a 
superposition of a finite number of certain selected modal 
functions. This type of analysis is often referred to as a 
Rayleigh type or the Rayleigh-Rit z method. Since a contin- 
uous system, such as a wing, possesses infinite degrees of 
freedom, the accuracy of the results depends, in general, on 
the number and choice of modal functions to be used. The 
two most common choices for these functions are either the 
coupled or uncoupled modes of oscillation of the system in 
a vacuum, since these functions satisfy structural boundary 
conditions. In determining the uncoupled modes, bending and 
torsion are assumed independent, such that the inertial 
coupling terms, those containing s a > are neglected in the 
equations of motion. Physically, this is equivalent to 
assuming, for pure bending, for example, that the mass 


(5-3) 

(5-4) 
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distribution acts along the elastic axis with no torsional 
deformation, and similarly for pure torsion, the wing is 
assumed to be constrained from bending. A coupled mode 
usually refers to a combination of torsional and bending 
deformations appropriate to the natural normal harmonic 
vibration of a freely oscillating system (30). 

On a wing with a true elastic axis and utilizing strip 
theory aerodynamics, it is computationally more convenient 
to use uncoupled modes in the analysis (20). Furthermore, 
Ref. 31, in comparing the coupled and uncoupled methods 
with experimental results on a straight, uniform, canti- 
lever wing with a large mass placed at various spanwise 
locations, indicates that using coupled modes analysis 
shows no better agreement with experiments, and that using 
uncoupled modes shows a more systematic improvement of 
accuracy when additional degrees of freedom are used. For 
these reasons, the uncoupled mode method will be used in 
this analysis . 

The number of modal functions assumed in any flutter 
analysis generally depends on the complexity of the struc- 
ture under consideration. Experience has shown, for simple 
wing models without large concentrated masses, that satis- 
factory results are obtained with two or three modes. The 
assumption is made in this analysis that the CMG mass is not 
large relative to the wing mass, otherwise little benefit in 
weight savings would result using active control. The CMG 
does, however, generate a concentrated torque on the wing, 
whose effect may require the addition of more modes in the 
analysis for sufficient accuracy. The procedure adopted in 
the following open-loop analysis will be to assume four un- 
coupled modes, consisting of the first and second bending 
and torsion modes of the uniform wing. This results in an 
18th order system. To determine whether the complexity of 
this formulation can be reduced without loss in accuracy, 
these results will be then compared with a two mode, 10th 
order, analysis. 

The bending and torsion deflections can now be written 
as 


My,t) = f hl (y)hi(t) + f h2 (y)h 2 (t) (5-5) 

a (y>t) = f a\ (y)c*i(t) + fa 2 (y) a 2( t ) * (5-6) 
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where fhl and fh2 are the first two uncoupled bending inodes, 
and f a i and f a 2 are the first two uncoupled torsion modes 
of a uniform contilever wing oscillating in a vacuum. The 
mode shapes are illustrated in Fig. 5-3 and can be repre- 
sented mathematically by 


fhl(y) » cosh(1.875*y)-cos(1.875*y)-.734|sinh(1.875*y) 

-sin(1.875.y)j (5-7) 

fh2 (y)=cosh(4 . 694 ’y) -cos (4 . 694*y)-l . 018^sinh(4.694*y) 

-sin(4.694*y)l (5-8) 


f «l(y) = sin -|y 
f a2(y ) = sin^- y 


(5-9) 

(5-10) 


To develop the ordinary differential equations of motion, 
the classic approach of Galerkin (32) will be used. Thus, 

Eq. 5-1 is multiplied by fhl and fh2 an d Eq. 5-2 is multi- 
plied by f a i and f a 2, and the resulting four equations are 
integrated over the length of the wing. Assuming two-dimen- 
sional strip theory, with the lift and moment per unit span 
being represented by the expressions defined in the typical 
section analysis of Chap. 2 , the equations of motion can be 
written, 

mFhlhi+^fhl (y G )hl+ F hl c hl h l + F hl,.alS*oi + ^ f hl(yG> f al(yG) a l + 

MGfhi fh? (yr:)h? .+ Fhl,a2 s a®2 + £®£fhl (yG) f a2 (yG)°2 

1 ? i r • i „ • i (5-H) 

= -77pb 2 |Fhlhi + U |F h i ja i a l + Fhl,a2 a 2j 

- ba^F hl>al a 1 + F hl >a2 a 2 ] } ~ 2 n P b u jj^ 

[ F hl b l +U(I' h l > al a l +F hl,o2 a 2 )+bft ' a ) (F hl.'Jl i l +F hl,a2a2)] 
x ♦(t -t)dt 

*NOTE: To be consistent with the typical section analysis, 

h and a are retained to represent the bending and torsion 
deflections. However, similar subscripts between the two 
analyses may refer to somewhat different quantities. Hope- 
fully, the formulation is straightforward enough so that no 
confusion results. 
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FIG. 5-3 UNCOUPLED BENDING AND TORSION MODES OF 
VIBRATION OF A UNIFORM CANTILEVER WING 
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+ -J^ f h2 f a 1 (YG ) ' “ l +s a F h2 , a 2 ® 2 + “]p f h2 f a 2 < y G 2 
= * npb 2 | F h2 h 2 +U ^F h2>a l a l +F h 2,a2 a 2] 

- ba [ F h2,al ci l +F h2 > a2 S 2jJ“ 2 77 P bu j 

[ F h2 b 2 +U (F h2 ,a I a l +F h2 ,a 2 a 2 )+b * a) < F h2 , al"l +F h2 , a2 a 2 >] * 

*(r-F)dF (5-12) 


s aGfr ^ ^ ^ , s aGi 


s a F hl 5 aA + -J^ f hl f al CyG) b l +s a F h2,alh 2 +_ f~ f h2 f al (yG) h 2 +I a F al a l 
2 

+ ^ f al (y G )5 l + ^ £ al f a2 (F G )5 2 +c al F <H a l + X f al ^G*" 

- "P b2 J ba j^ F hl,al’'l +F h2 > cil*'2j -U b (% _ a)| F ai®i" b2 0/«+ a2 ) F a l“l J 


+2np\h (ipha) 


y [ Fhi >° 


lhl+Fh2 , alh 2 +UFal a l+l > (%-a) F al^l * ( T -t )dt 


(5-13) 


5 aG, 


• » U/y P H *• j. . 

S a F hl,a2 h l + ~ f hl f a2 (yG) h l +s a F h2,a2 h 2 + “ f h2 f a2 (yG> h 2 +I a F a2“2 

+ ^5. £ a2yG^2 + -r^- f a2 f al( y G^l +c a2 F a2 a 2 + ^ £ a2(yGy 

l l t 

“ npb 2 j baj^Fhl } a2hl+Fh 2 ,a2h 2 j -Ub(%-a)^ a2 a 2 -b 2 (y 8 +a 2 )F a2 a2jJ 

+2jrpUb 2 (\ +a) j ^-^Fhl }a2 hi+F h 2 >a2 h 2 + UF a 2 a 2 +b (%" a ) F a 2 “ 2 j*( T ‘^) d t 

(5-14) 



- H w[ai f al<yG) +,i 2 f a2(yG)J +1 y° = T c 


(5-15) 


where 


F hl = 
F h2 = 
F a l = 
Fq2 - 
F hl, al * 
F h2 ,al = 
F hl ,a2 = 
F h2,a2 = 
C hl = 


u h2 

c al 

Ca2 


ii f hi 2 (y> d y 

j- r f h 2 2 (y) d y 

if f “ l2(y)dy 

n f ° 22(y)dy 

yf f hi f ai(y)dy 
i h 

j[ f h2 f ai<y) d y 

1 Jo 

jf f hl £ a2(y)dy 

1 cA 

H f h2 f a2(y> d y 

i Jo 

EIf hi f hi(y) d y 

if Elf h2 f h2(y) d y 

i h 


= _ i 

A) 
= _ i 

A 


GJ 


if 


£ ~ XX (y)dy 


al al 


GJ f a 2 f a5 I (y) d y 


(5-16) 


Another often used procedure for obtaining the equa- 
tions of motion is to first develop expressions for the kinetic 
and potential energies of the system and then apply Lagrange's 
equations (33). The resulting equations would, of course, 
be the same. 


5 . 3 Open-Loop Characteristics of Wing-Controller System 

Paralleling the procedure of Chap. 2 to determine the 
characteristic equation, Eqs . 5-11 and 5-12 are divided by 
7rpb^, Eqs. 5-13 and 5-14 are divided by 77 pb^, and the Lap- 
lace transform of the resulting equations is taken. The 
equations become, in matrix form, 
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|(> c allM- a )s 2 + l+(l*2a)*>(5)^ us 

+ 2 u 2 #( 5 )} F hl>a i 

+2ui(s) } F h2>a i 

{ jralM +(|+ a 2 ) ji 2 

+ a|%-a-(Jj-2a 2 )<i>(S)js 

+ r a 2 ^- u 2 ( l + 2 a ) o ( s )| F ot i 

2 -2 
r aG fia 


■ - Goi ® 


{( x al2M*a)s 2 + 1+ a-2a)v<s )Jus 
-»-2u 2 <f>(s)j Fhl ( c«2 

|( x «22A*- a ) s2 +^l+ (1— 2a^»<s)jui 
+ 2 u 2 *( i )| F h 2 f a 2 

-2 

W s 


{ [ r « 2 ^ + < l +a2) ] |2 

+uEs-a-(Ss-2a 2 )»>(i)j s 
+ r a / i .- u 2 ( U 2 aMs ) ^ F ^ 

“ g 02 s 



( 5 - 17 ) 



where, in this case 


"al 

u » TT-~- 

b oj 

*1 

In determining the effect of the CMG on the dynamic 
characteristics of the wing, the values of the nondimen - 
sional wing parameters corresponding to those previously 
used in Section 2.2_have been selected. The effect of 
the CMG parameter, G2, on flutter velocity and flutter 
frequency is illustrated in Figs. 5-4 and 5-5, respectively, 
for two different span locations of the CMG, corresponding 
to the midspan and wing tip. The CMG is assumed to be 
weightless. Also shown are comparisons between the two 
mode, four mode, and typical section analyses. 

It can be seen from these figures that the two mode 
analysis appears to provide sufficient accuracy for in- 
creases in the flutter velocity parameter up to and beyond 
3.0, and that the typical section analysis is seen to corres- 
pond closely to the three-dimensional case with the CMG at 
midspan. A comparison of the root loci, using the three 
methods of analysis, for particular designs having up = 2.5, 
is shown in Fig. 5-6. The effectiveness in increasing the 
flutter velocity by moving the CMG to the wing tip is evi- 
dent from Fig. 5-4 and is shown more explicitly in Fig. 5-7 
for various values of CMG weight. Increases in weight, 
generally decrease the flutter velocity, except in the 
vicinity of the wing tip, where the flutter velocity varies 
only slightly with CMG weight. Although the flutter velo- 
city remains fairly constant. Fig. 5-8 shows that the 
flutter frequency tends to decrease with CMG weight increases. 

From the above analysis, locating the CMG controller 
at the wing tip appears to provide the most effective control 
on the system. In Chapter VI, the feedback characteristics 
of various designs with this CMG location will be determined. 
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FIG. 5-5 VARIATION OF FLUTTER FREQUENCY WITH 
FLUTTER VELOCITY 


- 58 - 



IMAGINARY 



FIG. 5-6 ROOT LOCUS COMPARISON OF 3 METHODS 
OF ANALYSIS, up = 2.5 
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VI. FEEDBACK CONTROL OF THE THREE-DIMENSIONAL 
WING -CONTROLLER SYSTEM 


6.1 Introduction 


Obtaining solutions to the feedback control problem 
of a continuous, elastic, unstable system, such as a wing 
in an airstream travelling faster than the critical flutter 
velocity, can be most difficult. Exact solutions of just 
the aeroelastic problem have been developed for only a few 
special cases. Since the wing is actually an infinite- 
degree-of- freedom system, and thus has infinite states, 
new, more realistic, meanings must also be found for the 
concepts of controllability and observability (34). In 
Chap. 5, it was shown how the continuous system could be 
approximated by an equivalent system with finite degrees 
of freedom, and furthermore, in our particular case, how 
describing the system with only two modes of vibration 
can result in reasonable accuracy in determining the 
system’s flutter characteristics. There is no guarantee, 
however, that this simple model will be valid for the feed- 
back case. Higher modes may be excited and possibly driven 
unstable in attempting to control the lower modes. Thus, 
certain questions arise in examining the control problem 
of even the approximate system: 1) What is the minimum 

number of modes to be included in the analysis to insure 
that the controlled system is stable over the required 
range of velocities?, 2) What is the minimum number of 
controls needed for controllability and/or stability?, 

3) What are the essential feedback parameters to achieve 
stability?, and 4) How many measurements are needed to 
establish observability for the system? These questions 
will be discussed in some detail in the following sections. 

6 . 2 State Variable Selection and Controllability 

The method of design using arbitrary dynamics and the 
concept of controllability were discussed in Chap. 3. For 
convenience in determining the control characteristics of 
the three-dimensional system, the system equations, Eqs. 
5-17, are first put in standard state variable form. As 
in the typical section case, a decision as to which quan- 
tities are selected as states must be made, with the added 
complexity for the three-dimensional wing in deciding how 
many vibrational modes to assume in the analysis. Each 
additional assumed mode increases the number of states by 
four, due to the particular modeling of the incompressible 
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aerodynamics . Thus , for the two mode analysis, the system 
could be represented by ten states (two states corresponding 
to the CMG) , and eighteen states would be needed if four 
modes were included in the analysis. It is obvious that 
the amount of analysis and various design requirements such 
as sensor and feedback selection, could be substantially 
reduced if only two modes are necessary to obtain accurate 
control information. For that reason, and because the two 
mode analysis proved adequate in determining the open-loop 
system characteristics, a feedback control design using two 
modes was first attempted. The gains obtained from this 
analysis were then used to determine the characteristics 
of the system, approximated by four modes, to verify the 
accuracy of the two mode method. 

The selection of the ten states defining the two mode 
system is made consistent with those chosen in the typical 
section analysis of Chap. 3. They are 


h ll * h l 

(6-1) 

h 12 = h ”l 

(6-2) 

<*11 = a l 

(6-3) 

<*12 = “1 
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a 2 ~ b 

m 
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Following the procedure developed in Section 3.1. the 
appropriate portions of the set of Eqs . 5-17 are multiplied 
by D^(s) and the inverse Laplace transform taken, remem- 
bering that in the two mode case 

f h2 = f a2 =, F h2 * F a 2 = F h2,al = F h2,«2 = 

F rx2,hl “ F q; 2 ,h2 = ° 

The resulting equations are then written using the state 
variable definitions, Eqs. 6-1 to 6-10, and put in the 
matrix form, Fig. 6-1, where the matrix elements are 
defined in Appendix D. 

As discussed in Section 3.3, a necessary and suffi- 
cient condition for modal controllability is that the con- 
trollability matrix be nonsingular. It can be shown that 
this matrix defined by 

C = [g I FG : :F n_1 G 

is non-singular for the parameters selected in this analysis 
with Gn, Goi > 0 , u > 0. Furthermore, the system remains 
controllable, using just one controller, even when the num- 
ber of modes in the analysis is increased, so long as the 
quantities, an< * F ai^ y G) do not ec l ua F zero. This 

requirement is met with the CMG at the wing tip. These 
results are, of course, theoretical, with the assumed modes 
and aerodynamic model being approximations. Experimental 
verification of controllability of any design would be 
necessary to insure that the control or sensors are not 
located at the node of a particular mode of vibration. 

6 . 3 Design by Arbitrary Dynamics 

To investigate the feedback characteristics of the 
system defined in Eqs. 6-11, a particular CMG-controller , 
having a G 2 -value of 1.13 and CMG-to-wing mass and moment 
of inertia ratios of .05, was selected for analysis. For 
flutter velocity increases up to 50% (i.e., Up = 2.55), 
this system will have only two critical roots. Based on 
the typical section analysis in Sect. 4.3, this system 
will be the design having minimum power requirements, if 
the critical roots are placed on the real axis. Assuming 
a control law of the form 
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FIG. 6-1 TWO MODE EQUATIONS OF MOTION 
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and utilizing Eq. 3-26 results in feedback gains, varying 
with velocity, shown in Table 6-1. Since the three- 
dimensional wing is stable without the CMG for u<1.7, 
only the feedback characteristics for u > 1.7 were con- 
sidered of interest. 


TABLE 6-1 

SCHEDULED GAINS FOR FLUTTER -SUPPRESSION 
TO u - 2.55 


u 

*1 

k 2 

*3 

k 4 

K 5 

K 6 

*7 

K 8 

K 9 

*10 

1.9 

.006 

.009 

.057 

.014 

.006 

.080 

.0 

.001 

.001 

.005 

2.1 

.006 

.008 

.052 

.014 

.006 

.081 

0.0 

.001 

.002 

.005 

2.3 

.005 

.006 

.048 

.016 

.006 

.084 

0.0 

.001 

.001 

.005 

2.5 

.0045 

.006 

.043 

.017 

.007 

.086 

0.0 

.001 

.001 

.004 


As was discussed in the typical section analysis of 
Sect. 3.4, when there are only two critical roots the gains 
remain fairly constant with velocity. The relative sizes 
of the gains indicate that full state feedback may not be 
necessary. A number of reduced state feedback systems, 
using constant gains and meeting the stability criterion, 
are possible, three of which are shown in Table 6-2: 

TABLE 6-2 


CONSTANT REDUCED GAIN DESIGNS FOR u = 2.55 


K 3 

*1 



£l 

= .043 

K x - .0045 

K 3 = 

.043 

k 5 

= .020 

K 3 = .03 

K 4 = 

.017 

k 6 

= .086 

K 5 - .04 

k 5 = 

.020 
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A comparison of the root loci obtained for the sche- 
duled gains of Table 6-1 and the constant reduced gains of 
Table 6-2, case a), which requires the feedback of a , a , 
and a , is shown in Fig. 6-2. The final selection of any 
feedback scheme will depend on sensor and power requirements, 
with any variation from the gains in Table 6-1 resulting in 
increased power. 

The simplest control design would be to just feedback 
a . However, the large power increase associated with 
this design (cf. Fig. 4-5) would probably make it infeasible. 

6.4 Higher Mode Instability 

Although a number of designs have been developed for 
suppressing flutter in the system defined by Eqs . 6-11, 
there is no assurance, as yet, that the gains determined 
will stabilize this system with higher modes included in 
the analysis. An example of the occurrence of such a 
higher order instability is presented in Ref. 35 for a 
flexible launch vehicle. To guard against such occurrences 
in a design, usually more modes are included in the control 
analysis than are necessary to accurately determine the 
passive flutter characteristics. For example, in Ref. 36, 
although three modes were sufficient to predict the critical 
flutter velocity, seven modes were used in the control 
system design analysis. 

To determine the effect on higher modes of the designs 
developed in the previous section, constant gains, corres- 
ponding to u = 2.5 in Table 6-1, were used, and the second 
bending and second torsion modes were included in describing 
the wing-controller system. The resulting root locus is 
shown in Fig. 6-3. From the figure, the problem of higher 
order instability appears not to exist. Similar results are 
found for the gains in Table 6-2. This is not the case, 
however, if different root locations are selected for the 
critical roots. Figure 6-4 shows the locations of roots 
corresponding to different values of the imaginary part of 
the critical root at a velocity corresponding to u = 2.5. 
From Fig, 6-3, it can be seen that, although the roots 
associated with the first bending and torsion branches re- 
main relatively constant as the imaginary value of the 
critical root is increased, the higher order roots do vary. 
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FIG. 6-2 COMPARISON BETWEEN SCHEDULED GAINS 
AND REDUCED CONSTANT GAINS 
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with one root eventually becoming and staying unstable. 
Similar results will occur, if other prescribed root 
locations in the two mode analysis cause the gain values, 
and thus the power requirements to be large. The gains 
in Table 6-1 were selected by minimizing the value of power 
required for control upon encountering an impulsive gust 
disturbance. This criterion caused the minimum perturba- 
tion in the non-critical roots and resulted in low gain 
values. If other optimum criteria are established, for 
example, by defining certain weighting matrices in the 
quadratic synthesis analysis discussed in Sect. 4-1, pro- 
blems in higher mode instability may occur. 

6 . 5 System Observability 

As was discussed in Section 3.5, the sensor problem 
could be reduced if all the states could be estimated 
from a measurement of o . The system, Eqs. 6-11, will be 
observable if the observability matrix, defined in Eq, 3-27, 
is non-singular. It can be shown that the observability 
matrix is non-singular for the parameters selected in this 
analysis for Gqi, G]_]_, >0, u>0. The system does not re- 
main observable with a single measurement, however, if 
higher modes are included in the analysis, even if fhi(yc) 
and f a i(yG) do not equal zero. The relationship between 
the modes (cf.Eqs. 5-5 & 5-6) requires additional measure- 
ments as the number of selected modes increases, since only 
«(t) (or h(t)) and not «i(t) and «2 (t) can be determined 
from a single measurement. This illustrates another advan- 
tage to the design of Sect. 6-3, which required but two 
modes for a suitable design. 



VII. CONCLUSIONS & RECOMMENDATIONS 
7.1 Conclusions 


An analysis, utilizing modern control techniques for 
designing active flutter suppression systems with a con- 
trol moment gyro as the controller, has been presented. 

The following conclusions can be drawn: 

1. The concepts of modern control theory provide a 
well-organized, systematic, and mathematically 
powerful approach to the design of active flutter 
suppression systems. Using the method of arbi- 
trary dynamics, feedback gains, which produce 
stable systems having satisfactory response 
characteristics, can be easily obtained in a 
straightforward manner over a wide range of 
structural, geometric, and aerodynamic parameters. 
The solution for optimal gains based on a variety 
of cost functionals can be incorporated into the 
design not only for the single input/single output 
case considered here, but can be adapted to multi- 
input/multi-output systems. If full state feed- 
back systems are not feasible, these methods 
provide starting points in obtaining reduced feed- 
back and subop timal designs. 

2. The control moment gyro has certain attractive 
design features which warrant its consideration 
as an alternative to aerodynamic surfaces for ac- 
tively suppressing flutter modes. Controller 
designs characterized by simplicity, quick re- 
sponse, and low power requirements can be developed 

3. Although the flutter suppression analysis in this 
thesis was accomplished using a control moment gyro 
certain design results can be applied generally to 
any controller . 

(a) Controllability can be achieved using only 
one, properly placed, controller for signifi- 
cant increases in flutter velocity. 

(b) When using arbitrary dynamics, noncritical 
roots should remain constant. 
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(c) For introductory analyses, relative feed- 
back gain size is a reasonable design 
criterion for evaluating prospective designs. 

(d) Optimization criterion, such as minimizing a 
quadratic cost function, should be used 
cautiously, with the possibility of higher 
order instability kept in mind. 

(e) Two-dimensional analysis can provide valuable 
insight in developing practical design pro- 
cedures and determining basic system charac- 
teristics . 

(f) Certain aerodynamic surface parameters, such 
as flap width and/or length, can play a role 
similar to G 2 in selecting optimum designs. 

7 . 2 Recommendations 


In the field of active flutter suppression there are 
many possible areas of further research, many of which are 
discussed in Ref. 4, Recommendations evolving more directly 
from the results of this thesis are the following, some of 
which are presently being pursued at Stanford: 

1. Application of the control moment gyro concept as 
primary control device or in conjunction with 
aerodynamic surfaces, to an actual wing model, to 
determine realistic sizing and power requirements. 

2. Application of the design techniques, procedures, 
and tradeoff studies outlined in this thesis to 
an aerodynamic surface control, to compare power, 
performance, and complexity of design. 

3. Establishment of a practical design optimization 
criterion, thereby providing a more thorough 
understanding into the cause and prevention of 
higher order instability, 

4. Development of suitable mathematical modeling of 
the aerodynamics in subsonic, transonic, and super- 
sonic flow over complex elastic structures to allow 
the broad application of modern control technipes. 
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5. Establishment of sensor and estimator require 
ments to determine the required feedbacks for 
stable and optimum designs. 
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APPENDIX A 


AERODYNAMIC FORCES ON A TWO-DIMENSIONAL 
WING SECTION IN UNSTEADY INCOMPRESSIBLE FLOW 


Critical to the application of certain design techni- 
ques used in this thesis, such as the root-locus diagram and 
state variable formulation, is the particular modeling of 
the two-dimensional, transient aerodynamic forces using the 
Wagner indicial function. Although this concept is familiar 
to most aeroelasticians , it is felt that a brief description 
of its derivation is warranted. 

First, the aerodynamic forces in two-dimensional, in- 
compressible flow on a surface oscillating in simple harmonic 
motion are considered. Using potential flow theory, Theodorsen 
(22) obtained results for lift and moment in which the solu- 
tion was separated into two components, corresponding to 
circulatory and non-circulatory portions of the flow field. 

The resulting equations can be written 

L = /rpb^h+Ua-baaJ +2ffpUbC (k) £h-Ua+b (%-a)a:j (A-l) 

M * ?rpb 2 ^bah-Ub(%-a)a:-b 2 (|+a^)aj 

+ 27rpUb 2 (a+%)C(k)p+U«+b(%-a)aj (A-2) 


The non-circulatory components are contained in the 
first set of terms in Eqs . A-l and A-2, and are also valid 
for arbitrary motion as long as it is small. In the second 
set of terms, C(k) is, commonly known as the Theodorsen 
function and is evaluated by solving certain complex inte- 
grals resulting from the potential flow formulation. The 
Theodorsen function is complex and can be written after 

integration (2) „ v 

H X (k) 


C(k) - F(k) + i G (k) - Hl W (k)+ iHo W>W ■ (A ' 3) 
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where H]/ 2 ^ and H 0 ^ 2 ^ are Hankel functions of the second 
kind. A complex polar plot of C(k) is shown below. 



F(k) 


Fig. A-l COMPLEX POLAR PLOT OF C(k) 


These results have been used extensively for flutter 
prediction for many years, and can still provide a source 
of relatively reliable reference data for large aspect ratio 
wings even though they actually have supersonic flutter 
speeds (20) . 

Theodor sen's formulation for C(k) is invalid, however, 
for convergent oscillations in arbitrary motion. One 
method for solving this problem employs a superposition 
(Duhamel) integral of the linear results for simple harmonic 
motion. 

Corresponding to the steady state case, Wagner (23) 
determined the lift on an airfoil considering the growth of 
circulation about the airfoil as it starts impulsively from 
rest to a uniform velocity, U. The form of this lift is 


L = 2 npb Uw <j> (t) 

where w is the downwash velocity, t is a non-dimensional 
time, and (t) is the Wagner indicial function. Wagner’s 
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function was originally calculated by using a vortex-sheet 
method (23) . The method described here somewhat parallels 
the approach in Ref. 24. 

Because of linearity, it is possible to determine the 
lift due to circulation for any arbitrary motion of the air- 
foil by superposition. In an interval dt the downwash w(t) 
increases by the amount, dw(t) dt. When dt is small, this 

dt 

may be regarded as an impulsive increment and the lift can 
be determined by convolution: 

L = 2b77Up|w(o) 4> (o)+j" cfr(t-t ' ) dw ^ t ^ dt J (A-4) 

It is a convenient fact that in incompressible flow 
the circulatory component of force on an airfoil is deter- 
mined by the downwash velocity at the 3/4-chord (cf Eqs. 

A-l & A-2) which has the form 


w 3 /4c = ^ d* U & + b (^-a)or 


(A-5) 


Thus, only one Wagner function is required to describe the 
arbitrary motion. The Wagner function is evaluated by 
as suming 

ikf 

w(t) = w Q e 

in Eq. A-4 and comparing the results to Theodorsen's form for 
the oscillatory lift, Eq. A-l, which can be written 

L = 2 77 b p UC(k) w 0 e lk ^ (A-6) 

After some manipulation of terms, the Wagner function can be 
written as 

~ /-OO 

*(t) - ~ l sin kt dk (A-7) 

and its behavior as a function of t is illustrated in Fie 
A-2. 
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FIG A-2 WAGNER INDICIAL LIFT FUNCTION 

Although $(t) has a relatively simply graphical form, it 
is not mathematically expressible in terms of well-known 
functions. Therefore, several accurate approximations have 
been developed (20). For convenience in taking Laplace 
transforms, a particularly useful one for application in 
this thesis, is the form 

*(t) ~ 1 - .165e"* 0455t - . 335e~ (A-8) 
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Appendix B 

GYRO EQUATIONS OF MOTION 


In the twin-gyro controller, the equal but opposite 
precession of the gyro rotors can be obtained by using two 
torquing motors or by using one motor with proper gearing 
between the two gyros. The development of the CMG contri- 
bution to the typical-section equations of motion follows 
the derivation in Ref. 19, in which one motor is used. 

For this analysis, the following assumptions are made: 

(1) The CMG twin-gyro controller is rigidly 
attached to the wing. 

(2) The spin speed of each rotor is assumed equal 
and constant with respect to its respective 
gimbals . 

(3) The rotor spin axes are principal axes. 

(4) The moment of inertia of the typical section- 
controller combination about the elastic axis 
(y) is independent of the gyro gimbal angle 

a (Fig. B-l). 

(5) The typical section rotates only about the 
elastic axis. 



FIG. B-l TWIN -GYRO CONTROLLER ROTORS 
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The total angular momentum of the gyro controller is 


h TG = H G1 + h G2 + % 


(B-l) 


where 

Hi^g - total angular momentum of the twin -gyro 
controller 

H q1 = total angular momentum of gyro 1, including 
rotor and gimbal 

H G2 “ tota l angular momentum of gyro 2 

H-p = total angular momentum of twin-gyro 
controller "fixed" elements. 


From Fig. B-l 




cos o x + (hg-^ sincr 

T ■ A 

+ Ilz^z 



where 


(B-2) 


hg^ = angular momentum of gyro rotor 1 
I;Ly = moment of inertia of gyro 1 about y axis 
Il 2 = moment of inertia of gyro 1 about z axis. 
Similarly, 

**G2 = ~ h g2 coso ^ + ( h g2 sinCT + 1 2y“)9 ’ I 2z^’ (B-3) 

also 


h f = I F«y 


(B-4) 


where Ip is the moment of inertia of the fixed components 
about the y axis. Combining Eqs . B-l, B-2, B-3, and B-4 
and letting h gl = h g2 , I lz = I 2z yield 

^TG = sin a + Iya) y (B-5) 
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where 


H 

w 




x y - X ly + J 2y + ^ 

Equation B-5 illustrates the features of the twin-gyro 
controller in which the momentum change will occur about the 
y axis and crosscoupling terms cancel. 

To obtain the contribution of the controller to the 
typical-section equations of motion, Eq . B-5 is differentiated 
with respect to time in an inertial reference frame and ex- 
pressed in the wing reference frame by the theorem of Coriolis 
as 


i w _ 

h TG “ H TG + " x H TG = M w /G* ( b ~6) 

Here the superscript i indicates differentiation with re- 
spect to time in an inertial frame; the superscript w 
indicates time differentiation in the wing frame; to is 
the angular velocity of the wing frame with respect to the 
inertial frame which is this analysis equals ; and 
is the reaction moment of the wing on the controller. 

Equation B-6 can be written in scalar notation since 
the momentum and rotation are about the y axis. It becomes 

H TG = Iy“ + H w "C 08 " = %/G . 0-7) 

An equation involving the control torque and the gyro 
parameters can be obtained by differentiating Eqs. B-2 and 
B-3 and applying the theorem of Coriolis. This results in 

i-r = (Iiz^+hgl o sincr) Jt+ (hg^_ <7 cos a +I^ya)y 

H G 1 

+ (Ilz^-hgi CL cos <7 )z = Tq (B-8) 
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h G2 


(-I 2z <^+hg2 


a sin a 


)£+ (hg2 o cos cr 


+I 2 y«)y 


4- (-I 2 z ^+hg 2 a cos a )z = T 2 ' (B-9) 

where f and T 2 are the torques on gyro 1 and 2 respec- 
tively. 


Subtracting Eq. B-9 from Eq. B-8 and taking the, compo- 
nent in the Z direction gives 


I„ 0 - H « cos a 
Gz W 


- T lz - T 2 z 


(B-10) 


where 


^Gz - I lz + I 2z 


T 


lz 


T 2 z 


torques about z axis on gyros L and 2 
respectively. 


Since the gyros are geared together, 


T 


lz 


~T 


2z ‘ 


Therefore, let 


T lz " T 2z “ 2T lz " T c 

where T c is defined as the control torque supplied by the 
torque motor on the twin-gyro system about the z axis. 
Equation B-10 becomes 

t Gz° ■ H w“ COS(T = T c* (B-ll) 

With the assumption of small gimbal -angle oscillations about 
<7 = 0, Eqs . B-7 and B-ll can be linear lized by setting 
cos <7 « 1 . The linearized equations are 

!y “ + 17 • Ww/G 

-Hw “ + I Gz S = T c • 


(B-12) 

(B-13) 



Equations B-12 and B-13 are now used in conjunction with the 
equations of motion for the typical section to describe the 
motion of the wing-controller system and are applicable for 
use in the three-dimensional case under the assumptions of 
Chap. 5. 



APPENDIX C 


TYPICAL SECTION OPEN-LOOP CHARACTERISTIC EQUATION 

The characteristic equation described in Chapter II is 

s 2 |A 6 s 6+UA5 i 5 +(u 2 A 4 j+ 2 + ^4 3 G 2 ) s ^+u (u 2 A3 ]+A 3 2 +A 3 3 G 2 ) s ^ 

+ ^u^A2 ^+u 2 A2 2+A 2 3+ (u 2 A 24 +A 2 5 )G 2 j s 2 
+ u [u 2 Aii+A^2 + (^ 2 Al3+Ai4)G2] s 
+ u 2 (u 2 A 01 +A 02 +A 03 G 2 ) | = 0 

where 

A 6 = (I+^)(^+|+a 2 )-(x^-a) 2 

a 5 = d 1 A 6 +^(l+4)“^(x a +a)]+n 2 [2^(r^+a 2 )+(4x a a-^)/i-%] 

= d 1 A6+A5 D +n 2 A5 N 

A41 = d 0 A6+diA5 Dr +-niA5 N +n2[l-M( 1 + 2 a+2x Q ,)j 

a 42 * [(sj) 2(r a^ + t +a2)+r a'' (1+ ' i) ] 

A43 = (i+m) 

A31 = doAsnfnjjAsN+nifl-M 1 - 1 - 2 ^ 2 ^)] 

A32 = d l A 42 + (^|/ x ^‘ a ^ +n 2 ^ 2 r aM"(^> x ^’ 2a2) ] 

A 3 3 = 2n 2 +(l+4-)d^ 

a 21 = n G ^l-/i(l+2a+ 2x a )j 

*22 - d o A 42 +d l(s^) 4(%-a)+n 1 ^2r t ^-^^(J r 2a 2 )j 

* 
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A 2 4 = 2ni+(l+/0d o 


** ■ © ! . 

A 11 * d o(^|) >*ft- a > +n 0 [ 2r ^-(^) M%-2a 2 )j- ni ^ 

x M(l+2a) 

A 12 = djA^ 

A13 = 2n 0 

a 14 = d l A 25 

*01 - - n c(s|) ^< L+2a > 

A02 = d o A 23 
a 0 3 = d Q A 25 

F 71 = Y-(B^+A3 D )B4<%+a) 

A 6 

F 73 ’ X7<Bl+A5D> r <^ 

A 6 

F?4 " X“< B 1 +A 5D> 

f 76 - 7- G 1 (B 1 +A 5d ) 

A 6 



F 77 - '17 { n 2 A 5N +d lA 6 - jBj+Bj^^+a)]} 


F 81 - !t|u 2 [(%+a)(2 ni A 5D +B 1 d 1 +B 1 )-A 6 d 0 ]-B 2 B 4 (V+ a )} 
F 82 = -(2n 0 u 2 +B4d 1 )u 

f 83 = ~ [2n 1 u 2 A 5D +u 2 B 1 (d 1 +l)-B 2 B 4 ]-2A 6 n o a 4 l 
a 6 L 


f 84 = ^{^[ 2n l^ 2A 5D +u2B l( d l +1 ^" B 2 B 4], 


- A 6* 2 [ d o +B 3 n o +2n l] } 
f 86 = |2niu 2 A 5D +u 2 Bi(di+l)-B2B4| 

f 87 = _ ~*^ 2 { d o A 6 +n l A 5N~( d l +1 ^ [ B 2 +B l^% +a ^]} 

+B 4 [b 2 (%+ a )+r 2 /i+|+a 2 
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APPENDIX D 


ELEMENTS OF TWO MODE, OPEN-LOOP DYNAMICS MATRIX 


The open loop dynamics matrix, developed using two 
uncoupled modes of vibration, is written in Eq. 6-11. The 
elements of that matrix are defined as follows: 


F(l,2) 
F(2 , 7) 
F(2,9) 
F(3 ,4) 
F (4 , 7 ) 
F(4,9) 
F(5,6) 
F(6,4) 
F(7,7) 

F(7, 8) 
F(7,9) 

F(8,l) 

F(8, 2) 
F(8,3) 
F(8,4) 
F (8, 7) 


» 1.0 

= (r £yj. + | + a 2 )F al /A 

* -( x all^- a ) F hl,„l/ A 
= 1.0 

■ -< x «ll^ a)F hl,al /A 

= (1 +/ i l)Fhi/A 

- 1.0 

= _G ol 

= . j[(l+Mi)cli+l]F(2,7)F hl +[(x a . llM -a)d 1 +l 

+ B 3 n 2 ] F(2,9)F hljC(1 J u 
= 1.0 


=- [( x all/*-a) dl+l+B 3 n 2 ] F (4 , 9 )F hl >al + 



= - [d 0 + B 3 n Q + 2 ni ]ii 3 F h l,al 
* "j[“ 2 (( L+ f‘l) d 0 +2ri ^ f ^(3^ ] F hl F< ^ 2 ’ 7 ' ) 

+ [( x all^~ a ^ d o +d l +B 3 n l +1 ] G2F hl,al F (2,9)| 


I 
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F(8,9) - 

F(9,4) - 

F(9,7) = 

F (9 , 9) 

F.(9 ,10) = 

F ( 10 , 2 ) = 

F ( 10 , 3 ) - 

F(10 ,4) = 

F(10, 6) * 

F(10, 7) = 

F(10,9) 


( . 


j [( x all/J-- a ) d o +d l +B 3 n l +1 -J u Fhl >a l F ( 4 > 9 ) 

+ [ u 2 ^( 1+^1 ) d 0 + 2 n 1 ^ + 4 - j ^ ij 2 ] F h iF ( 4 , 7 ) 


“ G 0 x Gil 

| [( x o : ll / 1 ' a ) d l “( 1 + 2 a ) n 2 ] F h i sa iF ( 2 , 7 ) 

+ [(r 2 cx i^+§+a 2 )dl+^- -(%-a 2 )n 2 ]F a iF(2,9) 


u 


- j j^(r 2 ^.+i+a 2 )di+^- - (% _ 2a 2 ) n 2 j F a i(4 ,9) 
+ [(%ll/x-a)dl-(l+2a)n2]Fhl,alF(4,7)|G 


1.0 

(l+2a)n 0 u 3 F hl jQ;1 
-[r a 2 ^d 0 u 2 = (l+2a)n 0 u 4 ]F al 
-| u 3 jM d 0 >(^* 2 a 2 ) n 0 -( l + 2 a ) ni ]+ r a 2 ^ diujF al 

" G ol G ll d l n 


Gnd 0 u 2 

- j [( x all/ i -a)do-(l+2a)nij u 2 F hl)a lF(2, 7) 

+ [ u 2 r ( x « llM - a ) d 0 + ^ 3 di -(%- 2 a 2 )ni 

L *• 2 j 

-(l+2a 2 )n 2 ]+r a >] F a iF(2,9)| 

- j^u 2 j^x a nM-a)d 0 +^l di-(%-2a 2 )ni-(l+2a)n2j 

+ r a yjF«iF(4,9)+ [(x a iiM-a)d 0 -(l+2a)ni]x 
u 2 F h l , al F ( 4 , 7 )! 
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Where 


- (l+^l)(r a ^+g+a 2 )F hl F al -(x a ilF-a) 2 Fhi j «i 
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